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Abstract 

We analyze the axially symmetric pendent drop as it occurs in three different physical 
settings: Problem A, constant pressure, fixed circular opening (the siphon); Problem B, 
constant volume, fixed circular opening (the medicine dropper) ; Problem C, prescribed 
volume, constant angle of contact with a horizontal plate. As examples, the following 
results are verified. For Problem B we show that if the opening is small enough to 
support a stable pendent drop with a bulge, then as the exposed volume is increased, 
stable pendent drops with both a neck and a bulge will be formed. For Problem C we show 
that with increasing volume the profile curves for the family of stable pendent drops will 
develop an inflection point before instability arises. 

Introduction 


We first analyze the equilibrium and stability criteria for each of the problems. 
Problem A 


Here the drop is to protrude downward from a fixed circular opening of radius f, held 
at the level u = u where u is the vertical coordinate with positive direction upward 
and u - 0 is the zero pressure level of tne fluid. ( see Figure 1) . If X is the 
exposed body of the fluid and is the liquid-air interface with A(i>l ) its area, the 
potential energy of the configuration is 


) = aA(^2 ) + p g / zdV . (1) 

° ''x 

a is the surface tension of the liquid-air interface, p is the density of the fluid, and 
g is the gravitational constant. The condition for equilibrium is that the first variation 
of the potential energy 6E (.Q,N) = 0 for all normal perturbations N of Q which vanish 
on the boundary. The Euler Squations yield 


2H = -ku on k - p g/a . (2) 


H is the mean curvature of the surface measured so that it is positive at the drop tip. 
By a suitable scaling we may assume that k = 1. The condition for stability is that the 
second variation be positive for all non-trivial normal perturbations . 

^^E^(n,N) > 0 for all N ^ 0, but N = 0 on STt. (3) 


Problem B 


As in P'Toblem A the fixed circular opening of radius F lies in a horizontal plane, but 
now t.' exposed volume is prescribed (see Figure 1) . Now the condition for equilibrium is 
that the first variation of the energy 6 Eq(Q,N) = 0 for all perturbations N, vanishing 
on the boundary and for which the first variation of the volume is also zero. By the 
method of Lagrange multipliers we obtain 


6 (Eq+ XV) (f2,N) = 0 for some constant X (4) 

and all normal perturbations N vanishing on the boundary. This yields the condition 


2H = -ku + )* , where k= pg/a and X ■ is a constant, (5) 

By a vertical tramslation of coordinates we may take X* to be zero, reducing (5) to the 
condition (2) , and with the vertical coordinate of the opening at the level u = u. 
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( 6 ) 




The condition for stability is that 

6^(Eq+ XV) (n,N) > 0 

for all non-trivial normal perturbations N, which vanish on the boundary and for which the 
first variation of the volume is zero. 

Problem C 

The drop is now pendent from a homogeneous horizontal plate. The potential energy is now 

E(fi) = Eo(n) - aB|x| (7) 

where B is a physical constant and |z| is the area of contact of the liquid with the plate. 
Setting the first variation equal to zero for all volume preserving perturbations yields 
the conditions 

a) 2H = -ku + A for some constant X, k = pg/a (8) 

b) cosa =6 

Here a is the angle of contact of the liquid-air interface with the horizontal plate 
measured interior to the fluid. Again we may choose k = 1, and by a vertical translation 
of coordinates may set A = 0, with the horizontal plate at a level u = u. Clearly it is 
necessary that |3|< 1 so that 0 ^ a ^ tt . There are no possible pendent drops when a = it 
so that we may consider 0^ a < tt ( -1< 3 < 1) . As in Problem B the condition for stability 
is that the second variation 6^(E + AV)(r2,N) be positive for all non-trivial normal 
perturbations for which the first variation of volume is zero, (see Figure 1) , 

Problem A Problem B Problem C 
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Figure 1. The various drop configurations. 

Description of the Pr ofile curves 

Suitably normalized, the differential equation for the profile curve whose surface of 
revolution represents the liquid-air interface (satisfying (2) with k 1) is 


a) 

r* (s) = cos 

r(0) 

= 0 

(9) 

b) 

u* (s) = sin lij 

u(0) 



c) 

ijj*(s) = -(sin ^/r) -u 

1^(0) 

= 0. 



The solutions to this system have been carefully studied by many people. In particular, 

I should mention the work of D.V7, Thomson*^ , F.Bashforth and J.C. Adams^, and recently 
p. Concus and R. Finn^. 

There is a unique solution, {r (s ,< ) ,u (s , k) (s ,k) } , to t-he system satisfying the initial 
conditions, r(0,<) = 0, u(0,k) = u© = -2<, ii^(0,K) = 0, where k is the mean curvature at 
the drop tip. The solution exists for all s and all < and is analytic in both variables* 

We note that u - 0 gives a solution and that reflection of any solution about the r-axis 
yields another sclution. Drops with Uq < 0 represent pendent drops and the solutions 
corresponding to Uq > 0 represent "emerging" bubbles. We now list other important 
properties of the family of solutions. 

1. For "small" Uq <0 the solution can be expressed non-parametrically with u as a 

function of r over the entire positive r-axis, and u(r) 'v* uq where Jo(r) is the 

Bessel function of order zero. 

2. There is a value uj (uj ^ -2.5678-) such that the profile curve with drop tip at uj 

attains a simultaneous vertical tangent and inflection point at (r{,uj) where rj ^ .91 
and u? ^ -1.1 . For 0 < r ^ r? the curve is convex while for r greater than*^ r? 

the curve may again be expressed nOn-par3unetrically u = u(r). 
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u ■ u(r) 


3. For 
for all r. 


^8 


< 


4 . For Uq 
where 0 < r, 
decreases to**‘a 


«5 


< 

< 

valui 


0 the solutions may be expressed in non-parametric for*.., 
the profile curves attain a vertical tangent at a point (rx,u^) 


r? and u^ < u* • The curves form a bulge at this point and r 


^ 2 '. 




increasing functions of*^ u^ for u, 
approaches -«> . 


fo^s a nec 


at (r 2 ,U 2 ) where u, < 0, r 
< uj with limit r^ = 0 and limit u 


and 
1 = 


are 

as 


5. For Uq << ui the profile curves form a sequence of bulges and necks until it 
crosses the r-axis with r* (s) and u* (s) both positive from which point on the curves 
may be expressed non-parametrically u = u(r) out to r = «> . 


6. The first inflection point on a profile curve with tip at Uq <0 occurs at a 
point (f,d) where 0 < 0. f and d are monotonically increasing functions of u^ 
for Uq < 0. If Uq < u*, so that the profile curve has both a neck and a bulge, 
then the first inflection point lies between the first neck and first bulge. 



Figure 2. Possible drop configurations 
Analysis of stability 

Our method for determining the stable configurations for each of the problems proceeds 
as follows. Take a given profile curve {r (s ,u (s ,?) ,4^ (s , k) } and let (f,u) be a point 
on the curve, r = r(s,K> and H = u(S,k) . The profile curve from the drop tip to this 
point generates a possible pendent drop whose exposed volume V, can be calculated. 

V = volume of drop = Trr(ru + 2sin $) (10) 

The volume gives us a fourth function of the parameters s and k, V = V(s,k). For each 
of our three problems two of the four functions are prescribed. This generates a mapping 
from the (s,k) " plane (the parameter space) into a two-dimensional "control” space. The 
analysis of this map determines stability for each of the problems. 

Stability for problem A 

Here the appropriate map is A(s,<) where it is defined by 

A(s,k) = ( r(s,<), o,k)) . (11) 

The "control" space for this problem is the (r,u) “plane. It is easily checked that the 
derivative of A , DA(s,<) is invertible when s equals zero. Let O be the set of all 
points in the (s,<) -plane where the derivative is invertible. 

Definition. 0^ contained in 0 ie that component of O in the parameter space which 
contains the line^ s = 0. 

Theorem 1 . Every point (s,ic) in 0- determines a stedDle pendent drop for Problem A. 

( i . e . the drop generated by the profile curve (r (s,^) ,u (s,P) ) for 0 < s < s) . Any point 

outside determines an unstable pendent drop for Problem A. 

This theorem is essentially classical. 


288 









It follows that the "control set" A(Og) , is an open set in the (r,u) -plane. A point 
(r,u) determines a stable configuration for problem A if and only if it is a member of 

the set A(03). We now wish to describe the regions Og and A(Os). Since A(Os) is 

symmetric about the coordinate axes, we may restrict ourselves to the case r ^ 0, u 4 0 

Theorem ^ 2 . Choose r > 0 and consider the profile curve (r (s ,ic) ,u(s ) , s > 0 . 

There is a smallest positive value such that (s.k) is in Os for 0 < s < § while 
($,ic) is on the boundary of 05. On the interval 0 < s ^ S, we have r* (s) = cos 41 
positive so that 0 < < (ir/2) . 

Therefore the corresponding profile curve (r (s ,u (s , k) ) 0 < s < § may be expressed 
in nonparametric form u » f(r,?), for 0 < r < f where f = r(fi,i?) and H ® f(f,?). 

The point (f,d) lies on the boundary of A (Os) . It is the conjugate point to the drop 

tip along this profile curve. 

Since r* (s) is positive, we may use r as an independent variable rather than s. 

Points (r,u)^on the boundary of A (05) are determined by the condition that the 

derivative DA(r,<) be singular where A(r,tc) = (r, f(r,ic)) and u = f(r,K:) is the 
nonparametric representation of the curve. This occurs when f^(r,K) equals zero. (i.e. 
the point (f,u) is on the envelope , of the family of profile curves. 

Theorem 3 . The first envelope , of the fcunily of profile curves u - f (r,K)for 

< > 6 (uq < 0 ) and r positive, is the graph of a smooth (analytic) function u = e(r) 
for 0 < r 4 ao where olq is the first positive zero of the zero order Bessel function, 

• This fvinction has the following properties. 

limit e(r) = -<«> as r O'** , limit e(r) =0 as r aq 

The derivative e* (r) is positive on the interval 0 < r < a© with e* (ao) = 0 and 

limit e* (r) = +<» as r approaches zero. 

The entire envelope is thus a smooth curve without self-intersections which possesses 
a cusp only at (ao# 0 ). (see Figure 3 ) 

C onsequences. 

I. The mapping A(s,<) is a dif feomorphism of Os onto its image A (05) . 

II. For (r,u) in A (Os) where F < r* , the profile curve for the stcJDle pendent 
drop is convex. For r near a© ( F < Uq ) and (r,u) in A (05) , the profile curve 
will contain an inflection point, and so the stable pendent drop loses convexity. (Figure 3 ) 

III. There are no "inaccessible" stable pendent drops for Problem A. The vertical 
line r = r ( F less than ao ) intersects A (05) in a connected interval. Thus the 
stable pendent drop corresponding to the point (F,u) may be reached from the stcJDle 
zero pressure solution ( u = 0 ) corresponding to the point (r, 0 ) in A(Og) merely 
by increasing the pressure (p = -u) continuously from 0 to -u . (see Figure 3 .) 




Figure 3 . The envelope , stable configurations for problem A 
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stability for problem B. 


The constraints are now the rad as of the tube r, and the exposed volume V. Thus the 
control space is the (r,V) -plane and we are led to study the mapping B(s,k) from the 
parameter space to the control space defined by 

B(s,k) = ( r(s,K), V(s,k)) (12) 

As in problem A we let 0 be the open set in the (s,k)-plane where the derivative DB(s,k), 
is invertible. However, since B(0 ,k) = (0,0) the line s = 0 lies outside the set O. 

Fix P and consider the curve B(s,ic) for positive s. There exists a smallest positive 
value sg, such that the derivative DB(s,cr) is invertible for 0 < s < Sg but singular 
when s = Sg . Let (rg, Ug) be the corresponding point on the profile curve 
(r (s ,r) ,u (s ,^) ) , where rg = r(sg,ic) and Ug = u(sg,i?). It is a classical result that 
if (r,u) is a point on the profile curve prior to j^rg,u^) then the corresponding 
pendent drop generated by the profile curve up to (r,u) is "symmetrically" stable for 
problem B, while if the profile segment contains the point (rg,ug) then the generated 
drop is unstable for problem B. 

Definition. We call the point (rg,ug) the "Volume-constrained" conjugate point 
on the profile curve relative to the drop tip. 

Note. The axisionmetric pendent drop is said to be symmetrically stable if the second 
variation 6 "(Eq + AV)(vi,N) is positive for all non-trivial symmetric normal perturbations 
N, of Q which vanish on the boundary and for which the first variation of volume is zero. 
If the profile curve can be expressed non-parametrically in the form r = r(u) , then 
symmetric stability implies stability.^ In this case we observe that the angle of 
inclination si-, must be non-negat j '’e . However, if the angle of inclination becomes 
negative on some portion of the piofilc curve ( the corresponding drop is of re-entrant 
type) , then the drop is unstable for problem B due to a non-symmetric perturbation. 

This fact was noted by D.H. Michael and P.G. Williams For an alternative discussion 
see reference 5. 

Definition. Og is the subset of O consisting of all points (s,k) where 
0 < s < where Sg depends on It follows that (s,k) determines a stable _ 

configuration if it is in Og and an unstable configuration if it lies outside of Og. 

We now wish to describe Og and its image B(Os) contained in the "control space*, 
the (r,V) -plane. 

Theorem 4 . Og is a connected open set in the (s,k) - plane bounded on the left by the 
line s = 0 , and on the right by an analytic curve > , which is the graph of a positive 
analytic function Sg = a(^c). 

Therefore Og is c pen set in the parameter space and B(Os) is an open set in the 
control space. 

Theorem 5. Let (rg,ug) be the volume-constrained conjugate point on the profile curve 
(r (s ) ,u (s ,k ) ) . At the point derivative, r* (s) is positive. The point 

(rg,UB) is located between the first and second inflection points on the profile curve. 

If the profile curve possesses a bulge (and hence a neck) then (rg,ug) is located above 
the neck. As < approaches zero the point (rg,ug) approaches the point (a2^,0), where 
cii is a root of the equation rJo(r) + 2Jo*(r) = 0. (see Figure 4) 



Figure 4. Locations of the volume-constrained conjuoate point 
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By Theorem 4 the curve y is an analytic arc parameterized by k. Its image B(y) is a ! 

parameterized curve in the (r,V) -plane and, as in problem A, it is the en\^elope, 
the family of curves (r (s , k ) ,V ( s , k) ) . Thus Fg may be expressed in the form (rCK),V(K)) 
where r(K) = r(o(K),<) and V(k) = V(a(K),ic) are analytic functions of <- Furthermore 

limit (r(K:),V(K)) = (0,0) as <-► «, limit (r(K),V(K)) = (ai,0) as < 0 . 

By Theorem 5 we know that a given curve B(s,?), touches the envelope at a point where j 

(s) is positive. Thus, in a neighborhood of this point ^ express the J 

curves B(s,<) non-parametrically in the form V = g(r,K). if the envelope is smooth it j 

will be tangent to the family of curves V = g(r,K:), and itself would have a non-parametric j 

representation V = G(r). A point on the envelope of the family of curves V = g(r,K) , 1 

is determined by the condition g<(r,<) = 0 while the condition for smoothness is that \ 

gK:^(r,K) ^ 0. Since the envelope Fg = B(y) is a parameterized analytic curve it will i 

be smooth except perhaps at isolated points where the derivatives r* (k) and V’ (k) both 
vanish. At such points the possibility of a cusp arises. One such cusp occurs at (a]^,0). ! 

Conjecture. That part of the envelope Fg which lies in the half-space V > 0 is a j 

smooth curve which may be expressed non-parametrically in the form V= G(r) , 0 < r < j 

with G(0) = G(a]^) = 0 and G'(ai) = 0. There is a single value r* where G* (r*) = 0. | 



Figure 5. The curves B(s,k) and the^r envelope, Fg. 


If the envelope is a smooth curve, then it follows that the map B(s,K)is a diffeo- 
morphism of Og onto its image BCOg)* In this case any vertical line r « r in the 
control space would intersect BlOg) in a connected interval. The stable pendent drop 
corresponding to (^,^) is accessible from the flat drop u = 0 , corresponding to the 
point (F, 0 ) in the control space, through a smooth one-parameter family of stable 
pendent drops of increasing volume and fixed radius for the aperture until a maximum 
volume is reached. This is the proceedure used by E. Pitts in his paper. ^ 

If the envelope were not smooth and contained cusps, then the possibility arises that 
the mapping B(s,k) is not a one-to-one map of Og onto its image, or that for some f 
the intersection of the line r = f with the set : (Og) is not connected. In either 
case there would exist stable pendent drops corresponding to some control value (f,V) 
which could not be connected to (r, 0 ) in the manner described above. 

If we follow the "usual" proceedure of trying to describe those drops which are accessible 
from the initial drop u = 0 , corresponding to the point (r, 0 )then we have the theorem: 

Theorem 6 . (a) If r < rf then as we increase the volume from zero there will be 

produced a one-parameter family of stable pendent drops for Problem B. Through an initial 
range of volumes 0 < V < V-^{f) , the profile curves will be convex and the drops will 
develop a bulge. When the exposed volume reaches Vi(r), the profile curve will develop 
an inflection point at the edge of the drop. With increasing volumes the drops lose 
convexity, but before the limit of stability is reached pendent drops possessing both a 
necK and a bulge will evolve. 

(b) For f > where Ji( 6 j) = 0 the "drop" u 5 0 is unstable for problem B due to 
non-symmetric perturbations. For ? < the drop u s 0 is stable and with increasing 
volume the profile curves for the family of stable pendent drops will develop an inflection 
point before the limit of stability is reached. 

(c) For any radius f , drop height increases monotonically with volume throughout the 
range of stability. 

The result (a) of Theorem 6 was observed in the limiting case of small drop with 
narrow necks -bv A.K. Chesters 
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Figure 6. Drop formation for Problem B 
Stability for problem C. 

The control parameters are now the angle of inclination and the volume V, giving 

C(s,k) = (i|^(s,K) ,V(s,k)) (13) 

as our mapping from the parameter space to the control space* Again we let 0 be the set 
of all points (s,<)where the derivative DC(s,k:) , is invertible and we observe that the 
line s = 0 does not meet the set 0 since C(0 ,k) - (0,0). For each i? , there is a 
positive value Sp depending on such that DC(s,k) is invertible for 0 < s < Sp but is 

singular at Sp. Let Og be the set of all points in 0 of the form (s,^) where 

0 < s < s^. 

Theorem 7. 0^ is an open simply-connected set bounded on the left by the line s « 0 

and on tne right By a curve Yp which is the graph of a positive analytic function 
Sc = where limit is zero as k becomes infinite. 

Definition. For a given profile curve (r(s,i?)^u(s,?)) s > 0, the volume-constrained 
focal ^iht ^or problem C is the point curve with r^« r{Scr<) r u^«* u(s^,k). 

If a profile curve is to generate a physically meaningful configuration for problem C, 
it is necessary that the angle of inclination be non-negative along the segment of the 
profile curve generating the drop. Otherwise the drop would intersect the face. This 
eliminates from consideration the re-entrant drops. Therefore we let 0* be the set of 
those points (5,P) in Og for which the angle of inclination is positive^ 0 < s < 5. 

Theorem 8. The profile curve segment corresponding to any member of 0* generates 
a stable configuration tor problem C. If the point {s ,<) lies -outside the closure of 
the set 0* , then the generated drop is unstable for problem C. 

O 

Remark: This result which was essentially "classical" for problems A and B, is 

somewhat more difficult for problem C. 

In other words, let (rp,Up) be the volume-constrained focal point for problem C on 
some profile curve. Suppose that (r-,u^) comes before the point where the angle of 
inclination is zero. If (r.u) is a point on the profile curve prior to then 

the corresponding pendent drop is stadDle for problem C, while if (f ,Q) comes after 
then the resulting drop is unstable. 
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Theorem 9, The volume-constrained focal point (r for a given profile curve 

lies between the first and second inflection points, Tt comes ahead of the volume- 
constrained conjugate point for problem B, . 

Theorem 10. The set C(0*) is symmetric about the line ij; « 0. It is bounded on the 
leftTby tHe line ij; » 0. The^rest of C(0*) is bounded by where is the curve 

described in Theorem 7. ® 

As in problem the set C(Yp) is the envelope of the family of curves 

(i|>(s,k) ,V(s,ic) ) . By Theorem 9 each curve in the family will touch the envelope at a 
point where \j;*(8) is negative. Therefore in a neighborhood of the touching point 
each of these curves may be expressed non-parametrically V * h(4>,K:). The envelope 
is determined by the condition h * 0. It will be a smooth curve if h ¥ 0, 

If the angle is positive, theli dV/di|/ * h»^(4;,ic)= V /\l) is negative, and wfiere it is 
smooth, the envelope will be the graph of a decreasing®f unction. 

Conjecture. That part of the envelope , which lies in the first quadrant of the 
-plane is the graph of a smooth function V = V(\|/) 0 < < tt , with V* (0) * 0, 

V* (4^) negative for 0 < <-tt , and V* (tt) * 0. 

Computer calculations strongly support the conjecture. ^ If the conjecture is true, 
then (as in problem B) the map C would be a dif feomorphism of O' onto its image, and 
the intersection of a vertical line 4^ = ip , with C(O') would Be a connected interval. 
This would imply that as we move vertically along the line 4 = ip from (iT,0) to 
(if, V in the control space, the corresponding drops would generate the entire 

feunily *of stable pendent drops for problem C. If the conjecture were not true, then 
as was the case for problem B, the procedure just described would fail to pick up some 
stable pendent drops for Problem C. (see figure 7 .) 



Figure 7. The curves C(s,k) and their envelope 

The following theorem identifies those 8t^lble pendent drops that are accessible from 
drops of very small volume. 

Theorem 11. (a) For any angle of contact ip, 0 < ip < tt, stable drops of small volume 

are convex 7m d resemble spherical caps. These drops are generated by profile cu>“ve8 whose 
tip is at Uq, where Uq is large and negative. At a certain positive volume V, , where 
V. depends on ip^ the profile curve generating the drop will develop an inflection point 
at its edge. This drop is stable. As the volume is increased, further stable pendent 
drops are formed, and the inflection point on the profile curves will move to the interior. 
With increasing volume the limit of st^U^ility will be reached before a second inflection 
point appears. 

(b) If ip B 0 all profile curves corresponding to pendent drops of positive volume 
contain inflection point. Drops of small volume correspond to small values for u^ 

at tJie drop tip. As Up is decreased stable pendent drops of increasing volume are formed. 
The drop generated by that profile curve whose tip is at Uq * ^o * “2.5678- is unstable. 
(Computer results indicate that the steUt>le pendent drop of maximum volume occurs with 
^o * -1.6 with a volume of 18.4 ) 

(c) For any angle of contact drop height increases monotonically with volume 
throughout the range of sted>ility. 

Remark. For example, if the angle of contact 4 * (tr/2) , it follows that with 

bSlge®^?n instability is reached, pendent drops containing 
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Figure 8. Drop formation for problem C 
Justification of the stability criterion. 

As noted earlier, any re-entrant drop (drop for which the angle of inclination 
becomes negative) is necessarily unstable. Otherwise a stable symmetric pendent drop 
represents a strong local minimum of energy for any of the problems discussed. A n^ce 
proof is based of the method of H.A. Schwarz in his proof of the isoperimetric property 
of the sphere. For example, relative to Problem C we can show the following result. 

Theorem 12 . Let be a stable pendent drop for problem C with exposed volume V, _ 

and angle of contact 'i' , Consider any other pendent drop (Y,S) whose contained volume is V. 
For each horizontal plane p, below the supporting plane, let A{P) be the cross-sectional 
area of YOP, and let V(P) be the volume of the drop lying below the plane P. Let r(P) 
be the radius of the circle whose area is A(P) . Suppose that the pair (r(P),V(P)) 
determines a stable pendent drop for problem B and for every plane P. Then the energy 
E(n) (see (7)) is less than or equal to E(S), with equality only if (X,il) = (Y,S) . 

The method of proof is to first symmetrize (Y,S) producing a new drop of less energy 
and the same volume, which we then compare to the given stable drop. 
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